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Abstract 

Minijet production in hadronic and nuclear collisions through a 
BFKL pomeron ladder is studied for the energies of the future LHC 
heavy-ion collisions. We use unintegrated gluon densities compatible 
with the small- x increase of parton distributions observed at HERA. 
We show that at LHC energies the BFKL minijet and transverse en- 
ergy production is at most of the same order of magnitude as that in 
the collinear factorization approach. 
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Semihard parton scatterings with transverse momenta px ~ few GeV have 
been suggested to explain the rapid growth of the inelastic and total cross 
sections in high energy pp(p) collisions at y/s > 20 GeV |], |[. Especially in 
heavy-ion collisions at very high energies the semihard processes are expected 
to be abundant and dominate the transverse energy production in the central 
rapidity region [f3j. Event generators emphasizing the importance of semihard 
processes in ultrarelativistic heavy-ion collisions at y^i > 200 AGeV have also 
been actively developed during recent years [|], |]. 

Perturbative processes take place at the very early stages of the time 
evolution of an ultrarelativistic heavy- ion collision ||. In the central rapidity 
region, these semihard processes occur at time scales r ~ 1/pt, i-e. within 
the first fractions of fm/c, while the particle production from soft processes is 
expected to take longer, ~ 1 fm/c. Therefore, semihard particle production 
and associated transverse energy production give initial conditions for further 
space-time evolution of the formed partonic system, which eventually may 
lead to a thermalized quark-gluon plasma if the produced system is dense 
enough. The key feature of the semihard parton production is that it can be 
computed by using perturbative QCD || || . 

At central rapidities the semihard scatterings with momentum exchanges 
Pi ~ 2 GeV probe typically parton distributions at x ~ 2p T /y/s. Especially 
in nuclear collisions at the CERN Large Hadron Collider (LHC) with y/s = 
5.5 TeV, these x-values fall dominantly in the region of a rapid increase of 
the structure function F$(x, Q 2 ) as observed in deep inelastic ep collisions at 
HERA J?]. This increase persists down to Q 2 = 1.5 GeV 2 ||, thus strongly 
enhancing production of semihard partons as discussed in || . Nuclear effects 
in the parton distributions, in particular the observed nuclear shadowing of 
F 2 A (x, Q 2 ) at small x fpTOfl , are expected to be important, although the nuclear 



gluon distributions are not well known at the moment [|11], [12[ . 

In the above-mentioned studies, the production mechanism of semihard 
partons, minijets, is based on multiple independent 2^2 scatterings of 
partons. Collinear factorization is assumed to hold, enabling separation of 
the perturbatively calculable hard processes from the parton distributions 
containing non-perturbative input. In the leading twist approximation, the 
actual hard scattering involves only on-shell partons. In this paper, we will 
study a different mechanism for minijet production, which is not based on 
collinear factorization. 

The small- x increase of F^ measured at HERA can be explained by the 
Leading Log (logQ 2 ) DGLAP-approach |L3|], by the Double Leading Log 
(logQ 2 log 1/x) approximation [TJj], or by the Leading Log (logl/x) BFKL- 



approach [[Lq, [L6L 17|. In this study, we adopt the BFKL standpoint. Based on 



this we aim to study minijet and transverse energy production from a colour 
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singlet (cut) BFKL ladder spanned between the two colliding nuclei. We will 
not make any attempt to include nuclear effects here and, accordingly, our 
conclusions for the minijet cross sections should hold for pp collisions as well. 
The basic difference from the previous minijet studies including collinearly 
factorized 2 — > 2 on-shell parton scatterings only || , is that we now consider 
production of semihard gluons emitted by the virtual gluons in the legs of the 
BFKL ladder ordered strongly in rapidity, i. e. essentially from 2 — > 1 parton 
processes. This mechanism of minijet production was first suggested in 
T9fl . As a new feature, we use the unintegrated gluon distributions directly 



normalized to the integrated gluon distributions, which exhibit a small- x rise 
compatible with HERA results [20|j. Finally, we compare our results with 



the more conventional calculation M based on collinear factorization. Our 
conclusion is that the minijet and transverse energy production via the BFKL 
mechanism is at most of the same order of magnitude as the conventional 
one. 

Minijet production in hadronic collisions from a BFKL ladder in the case 
of two tagging (mini)jets separated by a wide rapidity gap was studied in 
21], UK and more recently in [|3|, |24|. A closely related subject in deep 



inelastic scattering (DIS) is the transverse energy flow due to minijets from 
the BFKL ladder |25|, p6| . Conceptually related to our study in this paper 
may also be the studies in [E7| , where collisions of two virtual gluon clouds are 
considered in the high energy heavy-ion collisions. Finally, let us also mention 



Ref. |p8fl , where minijet production by a usual soft pomeron is discussed. 

As a starting point for our problem, it is useful to recall inclusive minijet 
production in hadron collisions from a (cut) BFKL ladder as considered in 
[ [2"T|1 in the case of two tagging jets with a large rapidity gap in between. The 
situation on the parton level is illustrated in Fig. 1, where the transverse 
momenta of the tagging jets are denoted by k a x, kbT, and rapidities at the 
hadron CMS by y a , y b with y a ^> y^. The building blocks of such a ladder 
are the reggeized gluon propagators illustrated by the thick lines and the 
non-local gauge-invariant effective Lipatov vertices represented by the black 
blobs in the legs of the ladder. Squaring and summing all the gg — > (n + 2)g 
amplitudes in the leading log(s/i) approximation leads then in the multi- 
Regge kinematics to a hard pomeron, i.e. to a strongly rapidity-ordered 
perturbative colour singlet gluon ladder, denoted by /(q T ,k T ,y a — y b ) in 
Fig. 1. Note that in the multi-Regge kinematics, only the transverse momenta 
of the legs of the ladder become important, so that q 2 ~ — q T = — g T . Also, 
when intrinsic transverse momenta of the incoming partons are neglected, we 
have qx = — k a x and kp = k/r- An addition of a rung to the pomeron ladder 
is governed by an inhomogeneous BFKL equation O, jlBf. For details, we 



refer to the lecture notes of Del Duca [29|, here we merely cite the result for 
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the cross section, 



d 2 k aT d 2 k bT dy a dy b a ' ' d 2 k aT d 2 k bT dy a dy b 

where 

4iV 2 a 2 1 1 
d 2 k aT d 2 k bT dy a dy b N 2 - 1 fc^ T fc£r 

In the first approximation, only gluons can be considered, with densities 
g(x,n 2 ). The factorization scale is denoted by fi 2 . Fractional momenta of 
the incoming gluons are x a and x b , which in the multi-Regge kinematics 
become 

Vn k b T /0 N 

With the tagging jets the calculation is well defined, since large trans- 
verse momenta are required for the tagging jets so that the coupling of 
the ladder to the jets becomes perturbative. In this case one can still rely 
on collinear factorization and use the integrated parton densities xg(x,/i 2 ). 
Another important point in this calculation is that when all the radiative 
and virtual corrections are neglected, one is left with the inhomogeneous 
term of the BFKL equation only. At this limit the ladder shrinks into 
2/(qx, kx, y) — > 5^(qx — k-r), and the Born limit for the two jets separated 



by a large rapidity interval is recovered |29| . 

Within the framework of tagging jets, also more exclusive minijet pro- 
cesses have been studied by drawing a minijet out of the ladder p2[ . It is 



straightforward to show that after pulling out a minijet, the gluon emissions 
before and after the chosen minijet can be summed, resulting as a pomeron 
ladder on each side of the pulled minijet, as illustrated in Fig. 2. The cross 
section at the parton level becomes 



rfS _ 4:N 2 a 2 s a s N c 1 

d 2 k aT d 2 k bT d 2 k cT dy a dy b dy c N 2 - 1 tt , vr 



-2- f d 2 q 1T d 2 q 2T 



r(2) n . ^ 2/(k aT , qn, 2/a - I/c) 2/(q 2T , k iT , y c - ?/(,) 

(KcT _ QlT + Q2T7 To f2 ' VV 



h 2 h 2 



where the minijet drawn from the middle of the pomeron ladder is labelled 
as c. The factor a s N c / (7r 2 /c 2 T ) is a combination of the phase-space factor and 
the two Lipatov vertices associated with the step c. The factors 2 are related 
to the normalization of the ladder / as in Eq. (0). Again, if the radiative and 
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virtual corrections in the ladder are neglected, one recovers the Born limit 
for production of three jets well separated in rapidity |29[ . 



What we do in the following is that we simply relax the requirement 
of having tagging jets, since we want to study the leading BFKL minijet 
production mechanism, which is ~ a s . This means, unfortunately, that 
coupling of the pomeron ladder to the hadron becomes essentially soft in 
and &t, i-e. non-perturbative. This in turn results in the fact that while the 
other leg of the ladder is coupled to one constituent of the incoming hadron, 
in the squared graph the other leg may be coupled to another constituent, 



as shown in Fig. 3. This is nicely illustrated in a classic paper |TB[ where 
photon-photon scattering via heavy quarks was studied. Also, now that we 
do not require any tagging jets we have to give up collinear factorization, 
and, for the forward scattering amplitude, we do not have a perturbative 
Born limit to compare with, either. Therefore, the best we can do is to 



adopt the procedure for deep inelastic scattering in pOfl , where an addition 
of each rung into the pomeron ladder between the two hadrons or nuclei is 
expected to be described by the homogeneous BFKL equation. Justification 
for using a homogeneous evolution equation is that the ^-dependence of the 
hadron form factor describing coupling of the gluon ladder to the hadron is 
negligible with respect to the x-dependence generated by the BFKL evolution 

0- 

The homogeneous BFKL equation for the unintegrated gluon densities 
f(x,q^) given in [pOfl , which we will utilize in our problem, can be derived 



from the homogeneous BFKL equation for /(qT,kx,?/) in |T5| , [21]] by first 
integrating over the last momentum exchange kx in the pomeron ladder, 
then scaling with q\ and integrating over the azimuthal angle in the kernel, 
setting y = log(l/x) and finally scaling with x. The BFKL equation then 
becomes POfl 



df{x,q%) a s N c 2 /-°° dq( T 

~ x a = It / 

ox 7r Jo qf T 



2 r f(x,qj T ) -f(x,ql) f{x,q%) 



,2 „2 



+ 



Qt-Qit\ V9t + 4 9Tt 



• (5) 



The above equation is scale- invariant, so that additional information for 
fixing the normalization of the f(x, q^) is needed. The integrated gluon 
densities of a proton, determined through the global analysis of parton dis- 
tributions Pm, will provide us with this input: 



xg(x,Q 2 )= f ^f(x,q 2 T ). (6) 
j q T 

It should be noted that the non-perturbative nature of the pomeron-hadron 
coupling is now hidden into the initial (^"distribution f{xo, q^), which must 
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be supplied at fixed xq < 0.01. Equation (^) can then be solved numerically 
to obtain f(x,q^). In our discussion we will return to the treatment of the 
infrared problems of Eq. (||) but here will rely on the analysis of [^0| . 

Now we are ready to write down the formula for inclusive minijet produc- 
tion from the BFKL pomeron ladder between the two colliding hadrons (nu- 
clei) as illustrated in Fig. 3. Based on Eq. (f|), we get the factor a s N c /(7i 2 p\) 
when pulling out a minijet by fixing the momentum of one rung of the 
pomeron ladder. When doing this we form a new ladder on each side of 
the minijet. The coupling of these ladders to the hadrons or nuclei is con- 
tained in the initial condition for the unintegrated gluon distribution f(x, q 2 ), 
as explained above. Then, the cross section for inclusive minijet production 
is bound to have the same structure as in Eq. (f|), and one obtains (TTS], |T5[ 

d( ? ict j, ®sN c 1 /• 2 , 2 ( 2)f , s f(^i,qi) f(x 2 ,qi) 



K N ^-^\ /c? 2 qiT^ 2 q2T^ (2) (PT-qiT + q2T 

n 2 pi J 



d 2 p T dy " JV vr 2 f^i' 11 ' ' M " L ' qf T gf T 

(7) 

where pp and y are the transverse momentum and the rapidity (in the hadron 
CMS) of the minijet. From momentum conservation and multi-Regge kine- 
matics the momentum fractions become 

Xl ^^=e y , x 2 ^^e~ y . (8) 

Due to the fact that in this case we do not have an "external" hard probe like 
the virtual photon with an associated quark box as in DIS, nor an on-shell 
Born cross section to relax into, we cannot determine the overall dimension- 
less normalization constant Kn exactly. However, we are able to fix the slope 
of the minijet pp-distribution. Then, based on this slope we will discuss the 
upper bound for K N and actually argue that K N < 1. It is possible that an 
analysis of two-jet emission will clarify this issue. This demands in fact con- 
sidering a fcx-factorized form of the two-jet production cross section |3U|. By 
comparing this rate with the one evaluated directly from the BFKL ladder 
one would hope to get an absolute normalization to the off-shell gluon flux. 

Instead of numerically solving the homogeneous BFKL equation (|5D, we 
will use a simple parametrization for the solution, motivated by the form of 



the solution in the limit of asymptotically small x ||16| , p0| : 

a1la 2 \ 1 

(9) 



/( ^)=^(4v . a ex P r ln ^ 2 



xX ^(lV ^27rA" ln(l/ar) 



r2A"ln(l/x) 



In this expression A = 4a s ln2, A" = 28a s C(3), a s = 3a s /ii, ((3) = 1.202 is 
the Riemann zeta function and (p~Q, go an d q are parameters characterizing 
the initial distribution in [3(|. We reproduce the width of the solution with 
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an additional parameter r. For definiteness we shall take a s = 0.2, resulting 
in the same slope (~ x~ a5 ) as in the MRSD-' set of parton distributions [|31 



By choosing q — q — 1, C(p = 1.19 and r = 0.15, the proposed parametriza- 
tion (|D describes with a satisfactory accuracy the gluon distribution, which is 
obtained as the numerical solution of the BFKL equation in the global anal- 
ysis For comparison with Ref. p0 |, we show the unintegrated gluon 



distributions of Eq. @ in Fig. 4 as functions of q^ at different values of x. 

Proceeding then to the computation of the minijet cross sections, we inte- 
grate over the azimuthal angle between Pt and qi in Eq. (^) and fix = 1 
and y = 0. Results for the minijet cross sections da^ et / (dp^dy) at the LHC 
heavy-ion energy yfs = 5.5 TeV and at the UAl energy yfs = 900 GeV are 
shown in Figs. 5a and b. These are the main results of this paper. The cor- 
responding cross sections from the more conventional, collinearly factorized 
leading twist lowest order (CFLTLO) 2 — > 2 mechanism (see ||) with the 
MRSD-' parton distributions, are also shown for comparison. In the MRSD-' 
parton distribution functions we have made a scale choice Q = p^, so we can 
only proceed down to pr — V§ GeV, the lowest scale for this set. With the 
BFKL mechanism, we extend our computation down to 1 GeV. Note also 
that no .fT-factor to simulate the next-to-leading order (NLO) terms is used 
for the CFLTLO computation. 

In the BFKL approach only gluons are considered, and we have fixed 
a s = 0.2 in Eq. (^). In order to compare the two production mechanisms at 
the same level of approximation, we also plot the corresponding curves for 
the collinearly factorized results. However, since in the analysis of f(x,q^) 
in p0| , the strong coupling constant was set to run — phenomenologically, 



because NLO corrections to the BFKL ladder are not yet known — we show 
also the BFKL result with a s (pr)- 

Let us then come to the question of the normalization constant K^. We 
expect the BFKL mechanism to be valid and potentially important only at 
x < 0.01, where the rapid rise of the structure function F 2 (x, Q 2 ) is observed 
[0, ||]. Based on Eq. (||), the region of validity would then be < 9 GeV 
for -\/s = 900 GeV at y = in Fig. 5b. On the other hand, the experi- 
mentally measured jet cross sections at p T > 5.5 GeV can be explained 
through the collinear factorization mechanism for jet production rather than 
through the BFKL mechanism considered here. When the next-to-leading 
order corrections to CFLTLO jet cross section are included, the measured 
jet cross sections are well reproduced Let us also note that with a 

scale choice Q = Pt the NLO calculation can be reproduced by multiplying 
the Born level cross section by a factor between 2 and 1.5, depending on ^/s 

3 We are grateful to D.M. Ostrovsky for a numerical check of this statement. 
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and on the size of the jet |]34|| . This tells us that there is no room for an 
additional jet production mechanism of the same order of magnitude. Given 
this, and comparing the two mechanisms at the same level of approximation, 
the conclusion from Fig. 5b is that the BFKL mechanism in Eq. ([/]) can 
indeed at most have K N ~ 1. In fact, this provides an upper bound for the 
non-leading BFKL corrections as well. 

Up to this point, everything we have considered applies to pp collisions. 
The first estimate, without nuclear modifications to the parton densities, of 
the average number distributions of produced semihard gluons in central AA 
collisions can be obtained simply by multiplying the results in Fig. 5 by the 
nuclear overlap function, T^(b = 0) ~ A 2 /{nR A ) ||. For central Pb-Pb 
collisions, T PbPb (0) = 32/mb. 

As the last step, we estimate the transverse energy production due to the 
BFKL minijets in nuclear collisions. From the point of view of quark-gluon- 
plasma formation, and of its further evolution, we are mostly interested in 
the transverse energy deposit into the central rapidity region y ~ 0. An 
estimate for the rapidity distribution of the average transverse energy in an 
AA collision at an impact parameter b = is obtained as 

d£4 A m ,.,/". dai ct 



dy 



T AA {h) / dp T p T — — . (10) 
y =o Jpo dp T dy y= o 



Since the minijet cross sections computed for Fig. 5 are approximately con- 
stant at rapidities \y\ < 0.5, Eq. ( p~Q|) also gives a good estimate for the 
transverse energy deposit into the central unit of rapidity. By using the 
BFKL cross sections in Fig. 5 with po = 1 GeV we obtain the results shown 
in Table 1 for the LHC heavy-ion energy a/s = 5.5 ATeV. Again, comparison 
with the CFLTLO results with p = y/E GeV is made in this table. We 
remind the reader that the numbers quoted for the BFKL minijets are with 
Kj\t = 1 and we do not make any attempt to include nuclear modifications 



Li, 12] in any of the gluon distributions considered. 

Concerning our results, we would like to point out the following observa- 
tions: 

1. Although the BFKL mechanism we have considered here is leading 
in powers of a s — BFKL ~ a s and CFLTLO ~ a 2 s — and although the 
transverse momenta in the BFKL ladder have enlarged phase-space in the 
sense that they are not ordered as in the DGLAP-ladder, it is not enough to 
overcome the CFLTLO contribution at p^ > Po ~ 2 GeV. 

2. The slopes of the BFKL and CFLTLO computations in Fig. 5 are 
quite similar at large pt- This demonstrates that at fixed \fs the relevant 
physical scale is the transverse momentum p^. Since the gluon distributions 
have the slope ~ x -0 5 and since x ~ Pt/\/s, the jet cross sections scale as 
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dE? hPh /dy 

_L / ts 

y = 0,b = 


BFKL 
K N = 1, p = 1 GeV 


CFLTLO 
glue only, p = y/E GeV 


a s = 0.2 


3060 GeV 


3060 GeV 


« s = at s (p T ) 


4940 GeV 


4870 GeV 



Table 1: The average transverse energy at y — carried by minijets with 
Pt > Po in central Pb-Pb collisions at -\/s = 5.5 ATeV. 

d<ji et /dpT ~ (xiX2) _0 ' 5 /Pt ~ V^/Pt f° r both mechanisms. 

Let us next discuss the potential caveats and uncertainties. As discussed 
above, due to the non-perturbative coupling of the pomeron to the hadrons 
or nuclei, and lacking the corresponding partonic Born process as a special 
limiting case, we were able to fix the slope of the minijet production only, not 
the absolute magnitude. Based then on the measured jet cross sections []S2"| 
and the collinearly factorized calculations, we argued that the production 
of minijets through the BFKL mechanism can at most be compatible with 
minijet production in the collinear factorization approach ||. 

However, both approaches have problems when p^ — > 0, which we did 
not discuss above. At this limit the CFLTLO cross sections grow rapidly, 
and, especially with the HERA parton distributions, it is not plausible to go 
much below pt = Po ~ 2 GeV since this leads to an overprediction of the 
measured charged hadron distributions ||35|| . 

On the other hand, also the BFKL approach has serious infrared prob- 
lems. As seen in Eq. (|j|), a contribution to the BFKL evolution comes from 
the soft region < q^ ~ 1 GeV, the importance of which is enhanced if a 
running coupling is used. The infrared region was treated in [^(| by introduc- 
ing an x- independent form factor at q^ < q^ = 1 GeV together with a s (q\). 
For a more detailed discussion, we refer to [2(| and references therein. We 
believe that in spite of this phenomenological input in the soft sector, the 
correct order of magnitude of the BFKL results can be obtained since, first, 
the high-g!> region can be computed perturbatively from the BFKL equation 
and, second, the normalization of the unintegrated gluon distributions to the 
integrated ones can be made. We would also expect that if the BFKL mech- 
anism is indeed the reason for the rapid rise of the HERA parton densities at 
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small x, it could dominate the minijet production in the region of pt < Po, 
where the CFLTLO calculation does not apply. 

The subleading contributions are not yet known for the BFKL approach, 
although the work is currently in progress p6l. For the 2 — > 2 cross sections 



of massless partons, the NLO terms have already been known for a while [37 



and they have been applied to the case of jet production f33|| . However, one 



should keep in mind that minijet production is conceptually different from 
the jet production since minijets are not observed as individual jets. For 
minijets the NLO analysis can be made for infrared-safe quantities, such as 
the transverse energy deposit in \y\ < 0.5. This task has not been completed 
yet. Especially, it remains to be seen what the scale dependence of the 
transverse energy deposit turns out to be. 

The last potential caveat, well hidden in our BFKL approaclfl, is related 



to the inclusive and more exclusive processes. In pl| , the calculation is well 
defined, since the minijet production considered is fully inclusive. One finds 
a complete cancellation of infrared singularities between the virtual and real 
emissions in the pomeron ladder, leading to the explicit form of the kernel in 
Eq. (d). When the coherence of the ladder is broken by extracting a minijet, 
a potential problem with the exclusiveness and the cancellation of the sin- 
gularities arises (see f38|). At the moment we do not know precisely how to 
address this problem. However, since we are computing semi-inclusive quan- 
tities, in particular the transverse energy flow with less infrared-sensitivity, 
we believe that our estimates of the magnitude of minijet production with 
Pt > 1 GeV and of the associated Ejr- i n \u\ — 0-5 ar e reliable. 

Finally, we would like to discuss the connection to other studies. Recently, 
there has been a lot of activity in developing an approach to semihard par- 
ton production in ultrarelativistic heavy-ion collisions, based on an idea of 
colour charge distributions of the nuclear valence quarks moving along the 



light cone and creating clouds of virtual gluons around them [27]. In this 



approach, production of semihard quanta will result from a collision of two 
such incoming clouds. The leading mechanism of minijet production would 
then obviously be through 2^1 processes, i.e. of order a s , not a 2 s as in the 
collinearly factorized mechanism. Conceptually, this may be similar to the 
BFKL mechanism we have considered here, although, to our knowledge, an 
exact relation is not known at the moment. 

To summarize, we have computed minijet production from a BFKL pom- 
eron ladder and estimated the transverse energy production in the mid- 
rapidity in a central Pb-Pb collision at =5.5 ATeV. We have used un- 
integrated gluon densities similar to those of Ref. [|20||, so as to show that 



4 We thank Yu. Dokshitzer for pointing this out. 
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the BFKL contribution is at most of the same order of magnitude as the 
collinearly factorized one in the considered region. However, with the BFKL 
mechanism one could be able to bridge the way into the region px < Po in a 
manner at least partly controlled by perturbative QCD. 
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Figure Captions 



Fig. 1. Fully inclusive minijet production with tagging jets a and b |H 



The BFKL pomeron ladder is denoted by /(qx, kr, y a — yb), and only the cut 
ladder is shown. The thick lines represent the exponentiated gluon propa- 
gators and the black blobs stand for the effective non-local Lipatov vertices. 



Fig. 2. Production of minijet c from the cut BFKL pomeron ladder of Fig. 1 
[ p^j . A new ladder arises on each side of the minijet. 



Fig. 3. Production of minijets from the cut BFKL pomeron ladder without 
tagging jets. The ovals with one thick line and three thin ones attached rep- 
resent the incoming protons, and f(xi( 2 ), c^m)t) describe the BFKL pomeron 
ladders obeying Eq. For details, see the text. The cross section for 

minijet production through this mechanism is given in Eq. (|7|). 

Fig. 4. The unintegrated gluon distributions f(x, g T ) divided by ^Jq^ as 
given by the parametrization in Eq. (|D|) as a function of g T . The curves 
correspond to three different, fixed values of x. 

Fig. 5. a. Cross sections da^ et / (dp^dy) at y = as functions of for 
y^s = 5.5 TeV as predicted by Eq. (0). For the BFKL mechanism, two 
different curves are shown: the solid line corresponds to having fixed a s = 0.2 
in Eq. (|7|), the dashed line corresponds to a s = « s (pt)- For the collinearly 
factorized leading twist lowest order (CFLTLO) computation ||, we show 
three different curves: the dash-dotted line is the full computation with all 
parton flavours and a s (px), the dashed line is the line with gluons only, and 
the solid line is obtained with gluons only and a s = 0.2. To see the deviations 
from scaling limit (see the text), the log-log plot is used. 

b. The same as in Fig. 5a, but for y/s = 900 GeV. 
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